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Abstract 

To enable the study of criticality in multicomponent fluids, the standard spherical model is gen- 
eralized to describe an 5-species hard core lattice gas. On introducing S spherical constraints, the 
free energy may be expressed generally in terms of an 5 x 5 matrix describing the species inter- 
actions. For binary systems, thermodynamic properties have simple expressions, while all the pair 
correlation functions are combinations of just two eigenmodes. When only hard-core and short- 
range overall attractive interactions are present, a choice of variables relates the behavior to that of 
one-component systems. Criticality occurs on a locus terminating a coexistence surface; however, 
except at some special points, an unexpected "demagnetization effect" suppresses the normal diver- 
gence of susceptibilities at criticality and distorts two-phase coexistence. This effect, unphysical for 
fluids, arises from a general lack of symmetry and from the vectorial and multicomponent character 
of the spherical model. Its origin can be understood via a mean-field treatment of an XY spin 
system below criticality. 
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I. INTRODUCTION 



Criticality in liquid-vapor or fluid-fluid phase separation still warrants study: even after 
the advent of renormalization group theory, and its successful comparisons with experiment, 
open questions remain. One example is criticality in charged fluids such as electrolytes, 
molten salts, ionic solutions, etc. The long range of the Coulomb interactions impedes the 
application of most established methods and the interplay between electrostatic effects and 
long-range critical fluctuations is not fully understood theoretically. Indeed, the basic issue 
of the universality class of ionic fluids has been under debate for many years [l| and some 
questions still remain open. To gain insight into this and related problems, exactly soluble 
models can be valuable. Indeed, even if a model needs to be considered with circumspection 
in light of unavoidable simplifications, it may reveal significant features of criticality beyond 
those established by scaling and renormalization group analyses. 

In the history of models in statistical mechanics, the spherical or, equivalently, the mean 
spherical model [2, l3||, has played a special role. This "poor man's" Ising model has 
proved to be a mine of information because of its mathematical tractability: Thus only 

nriR 

as regards criticality, one can readily investigate ]4, 15|, 16| the role of dimensionality, scaling 
relations, finite size effects and the infiuence of long-range integrable interactions (such as 
l/r'^'^", where d is the dimension of the space and a>0). Consequently the spherical model 



has been applied in many physical situations, initially ferromagnets and 



quantum transitions 
forces 



j^, spin kinetics jl^, actively mode-locked lasers 



ater spin glasses 



critical Casimir 

12 1, etc. The model became all the more interesting when it appeared {3| that it 



belongs as a limiting case, n — oo, to the important class of spin systems in which n is 
the dimension of the order parameter (with n= 1, 2, 3, ■ ■ ■ for Ising, XY, Heisenberg, ■ ■ • 
models). 

It is natural, therefore, to consider spherical models with long-range Coulombic coupling. 
A pioneering investigation of a one-component plasma (OCP) spherical model has been 
undertaken by Smith but the limitations of an OCP model are well known and, in 
particular, a gas-liquid transition and corresponding critical behavior cannot be realized. 
Conversely, to treat electrolyte solutions a realistic model should first represent the neutral 
solvent, typically water; then two further species, namely, positive and negative ions, must 
be accounted for. Even if the solvent is appoximated by a uniform, structureless dielectic 
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medium, a colloidal system, for example, requires not only the macroions and their micro- 
scopic counterions but also the representation at some level of an ionic salt; thereby a ternary 
or quaternary system is called for. Accordingly it is desirable to develop spherical models 
for multicomponent systems. That is the aim of this paper. The investigation of the multi- 
component model proves interesting in itself although we will focus on the conclusions that 
can be drawn for simple binary fluids with short-range attractive interactions; applications 



to ionic fluids are presented elsewhere 



3, 
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Explicitly, we address a lattice gas with S species of particles, labelled a = 1,2, ■ ■ ■ ,S, 
in the grand canonical ensemble. Particles of a given species may occupy or leave vacant 
sites of only one sublattice so that the displacements separating the different interlaced 
sublattices introduce the crucial hard-core effects in a direct and transparent manner: see 
FigurelH We will use the vectors p = {pcr}, m, fi, h, etc., to denote the corresponding sets of 
densities, magnetizations, chemical potentials, magnetic flelds, etc., for the S species. Using 
the correspondence between lattice-gas and Ising spin models, and enforcing the S distinct 
spherical conditions with Lagrange multipliers, we extend to this multicomponent situation, 
the usual spherical model approach. This yields the free energy in terms of an iSx 5 matrix 
that describes the pairwise interactions between the different species : see Section [Til below 
and Eqns. (I22l) - (l28l) . It transpires that the singular part displays a form similar not only 
to one-component spherical models but also to Onsager's exact expression for the 2D Ising 
model [j, l5||. 

In the case of binary fluids (5 = 2), considered in Section UTTl major simplifications allow 
us to obtain explicitly all the thermodynamic and correlation properties in terms of the 
eigenvalues of the interaction matrix. The density correlations or (for charged fiuids) charge 
correlations for the different species, appear as combinations of two eigenmodes. These 
contributions become uncoupled only in the often considered but usually unrealistic fully 
symmetric case. 

To study critical behavior we go on in Section [TV] to consider systems with only hard 
cores and sufficiently attractive short-range interactions. With a proper choice of mixing 
coefficients, one can define m and m^, two linear combinations of the species mean magneti- 
zations or local densities, so that the usual spherical model critical singularities occur in the 
(T, m) plane at fixed Thence, a critical locus emerges in the full (T, fh, fh^) space which, 
together with a first-order surface, describes the infiuence of composition on the location 
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FIG. 1: (Color online) Illustration of a two-dimensional ternary lattice gas with species labelled 
0" = +, — ,0. Particles of each species occupy only one of the interlacing sublattices shown as dotted, 
broken, and solid lines, respectively; however, each particle may interact with all others via pairwize 
potentials cp^r^R'^ — R^)- 

of phase separation in the system. Via standard geometric arguments, criticality in the 
multicomponent model can then be deduced from corresponding one-component systems. 
Precisely, the same critical universality classes are realized as in short-range (attractive) 
spherical models. 

Nevertheless, a significant difference arises in the equation of state where a new, unex- 
pected mixing term appears. This can be understood heuristically as a type of demagne- 
tization effect arising as a consequence of the vectorial character of the model coupled to 
asymmetry and multispecies features. This term indeed suppresses the normal divergence 
of susceptibilities at criticality; furthermore, it induces a linear dependence of the chemical 
potential and pressure as functions of the total density in the two-phase region! These are 
certainly undesirable and unphysical features of any fluid model. This wayward behavior 
reinforces the remark [l^ that, because of the de facto vectorial character of the order 
parameter in spherical models, their predictions must be handled with perspicacity when 
modeling fluids. 

To gain some further insight into this unanticipated "demagnetization effect" we study 
in Section [V] an XY model beneath Tc using a mean-field approach in which the vectorial 
character of the order parameter, coupled to an asymmetry of the external fields, leads 
transparently to a very similar demagnetization effect. Finally, some general conclusions are 
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drawn in Section IVll 



II. MULTICOMPONENT SPHERICAL MODELS 
A. Fluid and spin systems 

We consider a (i-dimensional lattice fluid in the grand-canonical ensemble that consists of 
S species labeled cr = 1, ■ ■ ■ ,S, each being associated explicitly with only one of S identical 
interlaced sublattices (see Fig. [l|). Every sublattice is taken as the image of a periodic 
reference sublattice 7?.° after translation by a vector 5'^ so that every site i on a lattice a 
is characterized by a position = R° + <5'^. The reference sublattice is generated by the 
vectors aa (a = 1, ■ ■ ■ ,d), has a unit cell volume Vq and contains A/'= Hq, sites at positions 
R° = J2a -^i>^a specified by the integers -R°^ = 1, 2, ■ ■ ■ , A^^. 

It is well-known that a grand-canonical lattice fiuid is in correspondence with a canonical 
spin system Indeed, let us write the grand partition function of a fluid as 

S(T, M) = n ( E ^ E I - E ^'rNr)] , (1) 

where /3 is the inverse temperature l/k^T, while N^j and fi^ denote the number of particles 
and chemical potential of species a. The Hamiltonian H^'^ is expressed as a sum over 
particles = 1, ■ ■ ■ ,n'^ and / = 1, ■ ■ ■ , n"" as 

^'^^ = 1 E ^^rK-rJ), (2) 

(a,fc)7^{r,0 

where (/j^r is the pair interaction potential while is the position of the k-th particle 
of species a, occupying sites on the lattice a. Considering a system with hard cores, i.e., 
V5crcr(0) = +cxD, the sum in ([T]) refers to configurations where the local lattice density ncr(Rf ) = 
J2k ^(J^i ~ ^k) be only or 1, so that the local spin variable 

s.(Rn = 2n<,(Rn-l, (3) 

takes the values ibl as in the Ising model. A straightforward generalization of the procedure 
described in then leads to the partition function of a spin system, namely, 

S(T,/x) = n E exp(-/37^^P-), (4) 

(7 s,(R-) = ±l 



where the spin Hamiltonian is 

In this correspondence, the link between the coupling energies and pair interactions is 

J.,(Rr - RJ) = -\ ^..(R^ - RJ) if Rr 7^ RJ , (6) 
with Jcrcr(0) = 0, while the external fields are given by 

K = - l(pa, (7) 

with reference level 

0,= ^ ^.,(Rr„-RJ), (8) 
where is a fixed position. Finally, the background term in ([5]) is merely 

noi^^) = -l^^J2^''^ + l E ^^AR^-^])- (9) 

a {a,i)^{T,j) 

The correspondence between fiuid and spin systems follows straightforwardly for the other 
properties. For instance, the local density pa(j:)=ncr(Y) /vq is related to the local spin via 

v,p^{'R!l) = \MK1) + l] , (10) 

while the species correlation functions 

a,(Rr - RJ;T, A*) = (p.(Rnp.(RI)) - (p.(Rr))(p.(RJ)) , (11) 
are related to spin correlations via 

t;o^a.(Rr -RJ;T,Ai) = \ [(s.(Rn^5.(RJ)) - (s.(Rn)(s.(RJ))] . (12) 

As usual, the angular brackets denote grand-canonical expectation values. 

In order to define density or charge correlations simply in this lattice geometry, it is 
convenient to work in Fourier space. We consider periodic boundary conditions and define 
Fourier series with respect to the reference sublattice IZ^ by 

s.(k) = 5^e-^-^^s.(Rn. (13) 
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Then, when J^r is periodic over the reference sublattice, we may write 

X.(k) = Y: ^'-■^^^J.riB^ - RJo) = E e-^'^-(^--^I) J^,(R^, - RJ) , (14) 
* j 

with any fixed positions zq and jo- The wave vectors should be combinations of the reciprocal 

vectors (defined by aa-bo,/ = 27r5a^o') such as ^ = Y^a ^a^a with /cq, = 0, il/A^^^, ±2/Na, ■ ■ ■ . 

In the following, the first Brillouin zone is denoted as B. The density correlation function 

Givjv, and, for fiuids of particles carrying charges q^, charge Gzz and charge-density G^z 

correlations can then be defined via 

G,,(k; T,p) = J2 g^g^G..(k; T, p) , (15) 

cr,r 

where X and Y stand either for N or Z, with -^j^ = 0, i)z = ^- We also define structure 
factors as 

5,,(k; T, p) = -—^G,,(k; T, p) , (16) 

where q is an elementary charge, while the total density is p = YliJ\Po)- The term Vq 
compensates here for the homogeneity difference between the discrete and continuum Fourier 
transforms. 



B. Mean spherical model 

We are not able to perform the multiple sums in dH) in general. Instead, we adopt the 
appropriate mean-spherical model [3|| and compute the multiple integral 

S'(T,h)= jWds^e-^'^^ (17) 

with 

H' = H^P''^ + J]A.J]4(Rn- (18) 

a i 

As usual, the Lagrange multipliers X^j are introduced to allow imposition of the mean spher- 
ical conditions which need to be enforced uniformally for every species; specifically, the 
relations 

{Y.'l^^^))=^ ' ^ = l,---5. (19) 
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define the Lagrange multipliers or spherical fields as implicit functions of (T, h). Conse- 
quently, the free energy per site (of the reference sublattice) is 

-/?/[T,h,A(T,h)] =lnS'(T,h)/Ar, (20) 

in terms of which the spherical conditions (fTOl l can be rewritten as 

df 



(4) 



= 1 , a = l,...S. (21) 



As already remarked in SeclH the spherical model in this form describes exactly spin models 



with fixed 
scalings 
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ength or continuous n-component spins in the limit n — > oo with appropriate 

As standard Pl, the calculation of S' is performed in Fourier space. For consistency, we 
will suppose that J^r satisfies the symmetry condition Jo-T-(Rf — RJ) = Jar ~ I^J)] ■ 

The calculation is then a straightforward generalization of the mean-spherical techniques 
used for single-species systems. The free energy per site can be decomposed into a sum of 
three parts: f = fs + fh + fo- The singular part of the free energy is 

- PUT, h) = -^Yl 1^ {/^^Det [A(k; A)] } . (22) 
where the sum runs over the reference Brillouin zone B while A(k; A) is the S xS interaction 



matrix with elements 

A,,(k; A) = S^^r [X. + ^ZM] - 1(1 - 5.,r)ZAk) , (23) 
in which for any function ^(k) we employ the notation 

Ag{k) = \ [m - m] , (24) 
while, dropping the tildes in (fTSll and (|2T]1 . the shifted or net spherical fields are 

A. = A, - \ZM . (25) 
On the other hand, the /i-dependent part of the free energy is given by 

-/3/,(T,h) = i/5(h|A-i(0,A)|h), (26) 
while the analytic background part, following from Q, is 

- /3/o(T, h) = i5 In TT - l3nolM , (27) 



which will be neglected henceforth. Because of the logarithm in (1221) . these results are valid 
while the eigenvalues of the matrices A(k; A) are positive for every k; when one vanishes, 
the expressions f l22] ) and f l26l) become singular and phase transitions are implicated. 

The last step is taking the thermodynamic limit A/" ^ oo (valid provided the Fourier 
transforms remain well-defined) with the result 



PfsiT,h) 



ln{f3^Det [A(k; A)]} 



(28) 



where /j^ is a short-hand notation for j^^^ Vod'^k/ (27r)'^, while the h-dependent free energy fh 
is still given by (l26l) . At this point, it is worth noting that the structure of fs, as an integral 
over the Brillouin zone of the logarithm of the interactions in Fourier space, is similar to that 



present in Onsager's exact solution of the 2D Ising model 
systems are dramatic [15|, 116|. 



2l| . The consequences for charged 



17l | since this form determines the coupling or decoupling of 



correlations in symmetric and asymmetric systems. 

With these results in hand, we find that the mean particle densities pa = (pcr(Rf)) are 
related to the mean magnetizations via 



2p„VQ - 1 



is.) 



(29) 



which, in turn, enter the free energy in standard manner as 

fh = ■ 

a 

As a result of (l26l l and ( l29l ). the link between h and m is then merely 



(30) 



h = - i0 = 2A(0)m, 

where ^ = {pa-} and we have recalled (I7|) and introduced a fixed vector (t) = {(pa-}'- see 
Finally, in the thermodynamic limit, the spin-spin correlation functions are given by 

1 dU - Q 



(31) 



(32) 



T,h,A 



where 2 — d^^r is merely a symmetry factor, while we recall (fT2l) for the density correlation 
functions Gr„-- 
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III. BINARY SYSTEMS 



The previous analysis holds for an arbitrary number of species. From here on, however, 
we focus on the simplest case, i.e., binary mixtures with species labels 1 and 2. For many 
properties, it is useful to decompose densities, chemical potentials, etc., in terms of means 
and differences; so for every function (or g^a) we define 

9 = 1(91 + 92), 9^ = 1 (91 - 92) ■ (33) 

Moreover, for simplicity, we suppose that the translation vectors 5"^ = 'Yl,a^a^a satisfy 5a = 
or 1/2 so that the Fourier transforms J^t are real. 



A. Basic features 

For the case 5 = 2, simplifications allow more explicit results. First, let us introduce the 
energy scale 

Jo = 1^12(0) , (34) 

and, following ( 1331) , write 

A J(k) = \ (AJn + AJ22) , (35a) 
Ajt(k) = l(AJn-AJ22). (35b) 

Then the eigenvalues of the 2x2 matrix A may be written 

A±(k; A) = A + Aj(k) ±Z}(k;A), (36) 

where 



D(k; A) = ^ [At + AJt(k)]^ + i[Ji2(k)]2 > 0. (37) 

As remarked above, these expressions are valid when A_ and A+ are nonnegative while 
singularities arise only when A_(k; A) [< A+(k; A)] 0. 

Now the argument of the free energy integral in (l28l ) is ln{/3^Det[A]}, where the deter- 
minant of the interaction matrix can now be written 

A_A+(k; \)=u + 2AAJ(k) - 2A^Ajt(k) + Aj2(k), (38) 
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where we have introduced the crucial parameter 

m(A) = A-A+(0; A) = P - At' - jo' , (39) 

which vanishes when Det[A(k)] vanishes at k = 0, while the squared interaction term in 
fl38l) . namely, 

AJ^(k) = Jo' - i^il(k) + AJn(k)AJ22(k), (40) 

vanishes as |k^|. 

In terms of the eigenvalues, the spherical conditions ( l2Tll become 
Finally, the /i-dependent part of the free energy entails 

^-<-'^^ (;::)• 

while the magnetization-field or density-chemical potential relation (ISTll becomes 

= m^X„ - rririo , t ^ a. (43) 

At this point, Eqs. (i34l) -(l43l) entirely define the system and the need is to analyze their 
structure and consequences. 

B. Correlation functions 

The density pair correlation functions are given generally by (fT2l) and (l32l) which, when 
5 = 2, reduce to 

In terms of these, one can use ( fTSi l to obtain the overall density-density correlation function, 
Gjvjv and, the complementary compositional correlations or, for charged systems, the charge- 
charge correlation function Gzz- 
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From a purely mathematical perspective, it is also instructive to decompose the fluctu- 
ations with respect to the eigenvectors of A which, of course, depend on the wavevector k 
and the fields A. Thus if we define 0(k) via 

tan</)(k) = 2 {D{k; A) - [A^ + AJ^(k)] } /Mk), (46) 

it can be interpreted as the angle determined by the eigenvector associated with A+(k; A) 
relative to the a = l axis. Then if we introduce the density fiuctuations p"*"(k; A) and p^(k; A) 
via 

= [pi cos 0(k; A) T P2 sin 0(k; A)] / ^2, (47) 
and define the corresponding correlation function G±± in the natural way, we find 

while G'+_ = G-+ vanishes identically. 

However, the eigenmodes ( l47l) will rarely be of direct physical significance. Rather the 
physically accessible fiuctuations, represented in particular by the structure functions in- 
troduced in (|T5l) and (flGl) . will typically involve a mixture of the underlying eigenmodes. 
Specifically we find 

SMM{k-T,p) ^ B{k;X) 1 - B{k; A) 
k,T/Apvo A-(k;A)"^ A+(k;A) ' ^ 
and, for charged systems with q+ = —q^ = q, 

S,,{k;T,p) 5(k;A) , 1 - 5(k; A) 



KT/Apvo A+(k;A) A_(k;A) 
where the mixing amplitude B is 



(49b) 



i?(k;A) = i + iJi2(k)/D(k;A). (50) 

Evidently, singular behavior, anticipated at criticality in A„(k; A), will in general affect both 
zz as we discuss in detail elsewhere |l5l. Il6l. Il7l|. 
However, a special situation arises when the two species 1 and 2 are symmetrically related 
so that Jii(k) = J22(k) which implies, via (135bl) . AJ^(k) = 0. For a charged system this 
corresponds to complete charge symmetry as exemplified most simply in the restricted prim- 
itive model (RPM) of equisized hard spheres with charges of equal magnitude but opposite 
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sign. But neutral systems where species 1 and 2 differ only in chirality demand a symmetric 
description quite naturally. Then, on the locus of symmetry where pi = p2 (corresponding to 
electroneutrality in 1:1 ionic fluids) one has fii=fi2 and, hence, via (l3T! ). Ai = A2 and thence, 
via ( l33ll = 0. In this case one sees from ( l50ll that B{k; A) vanishes identically so that 
the eigenmodes precisely specify S^n and Szz which, therefore, become totally decoupled\ 
This turns out to play a crucial role in the study of charge screening near ionic criticality 



3^ 
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17| albeit for generally unrealistic charge-symmetric systems. 



A small technical detail deserves mentioning in this fully symmetric case if Ji2(k) should 
change sign for k 7^0 (which is not unreasonable); then the ratio Ju/DCk; A) in (iSOll together 
with A+ and A_ involve nonanalytic absolute values but in such a way that the combinations 
Snn and Szz in ( l49l ) remain completely analytic. 

Finally, the cross charge-density structure function is also expressible as a combination 
of the two eigenmodes via 

S^z{k;T,p) At + AJt(k) 



A_(k;A) A+(k;A) 



(51) 



k^T/8voP D{k;\) 

As is to be anticipated, this vanishes identically on the symmetry locus when (1, 2) symmetry 
is present. 



C. Appropriately mixed thermodynamic variables 

Depending on the symmetry of the system, the previous relations may be handled more 
or less conveniently. In the general asymmetric case (Jn 7^ J22), the spherical constraints 
f l4Tl) can be rewritten as 



and 



while the external fields are given by 



h =2 [m(A - Jo) + m^^'^] , (54a) 
=2 [mX^ + m\\+ jo)] • (54b) 
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Note that, since A_ and A_|_ are nonnegative, the condition (l52ll is consistent with the 
expectations |m|<l and |m^|<l. 

For further analysis it is convenient to introduce the basic integral functions 

which are simple generalizations of the typical integrals involved in the analysis of standard 
spherical models. 

Now, in the fully symmetric case, AJ^, A^, and vanish identically so that the rela- 
tions (l53ll and (I54b[l have no role to play. Then (l52l) and (I54ap closely resemble the basic 
expressions for the single-species (5 = 1) or standard spherical model. These in turn lead to 
the basic equation of state which, in terms of the reduced temperature variable 

t = (T-T,)/T„ (57) 

can be written most transparently near the critical point (T = Tc, m = 0) as 

Pou^^'' ^ Ctt + CraVr?, (58) 

where, recalling (l39l) . m = A^ — Jq while 7>1 is the fundamental dimensionality-dependent 
exponent, and po, Q, and Cm are fixed positive coefficients. 

Now physicochemical insight into the behavior of binary fiuid mixtures suggests strongly 
that their critical behavior will, when expressed in terms of suitable density and field vari- 
ables be essentially the same as for a single-component fiuid. However, the "suitable" or 
"appropriate" variables will, in leading order, be linear combinations or mixtures of the re- 
lated binary thermodynamic variables, specifically, the fields and densities. Furthermore, 
the appropriate mixing coefficients must, in general, be nontrivial functions of the state 
variables. 

It follows that our primary task now is to find what the appropriate mixing coefficients 
are. To that end, we introduce the general linear combinations 

A = I (^-^Ai + 6X2) , At = 1 {e-'X, - 6X2) , (59) 

together with corresponding remixed interactions 

AJ, Ajt = 1 {9-'AJu ± ^A J22) . (60) 
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The mixing parameter 9 is to be determined later. In terms of these new variables and 
interactions, the basic determinant becomes 

Det[A] =u + 2AAJ(k) - 2X^Aj\k) + AJ\k), (61) 

where, following (l39ll . the value at k = is now 

w(Ai,A2)=A2-At'-jo'. (62) 

Then it proves necessary to introduce a second state-dependent mixing parameter 6m by 
writing 

m,rh) = \ (6'~^mi ± d„,m2) • (63) 

In terms of these new variables and the integrals ( l55| l and (l56ll . the original spherical con- 
ditions fl4TI) become 

1 = e-^ksT{x - x^)gix) + C2ix)KT 

+6'^(m^+m^ + 2771771^), (64a) 

1 = ek^Tix + x^)gix) + jCii\)ksT 

+em{m^ + rh)"^ - 277i7v)). (64b) 
Finally, it is helpful to define new external fields via 

h = i + 9hh,) , /it = i (0, 1/^1 - 9f,h2) , (65) 

which are linked to the generalized magnetizations via 

h = 77l[{X- jo) UJ+ + X^iU- + joTT+j 

+ mt [(A - Jo) ^- + At^+ + jon+,] , (66) 

and similarly for h'^ in terms of 7r~ and Tr"^ with and interchanged while the coefficients 
iu±, TT^, and TT^t are found to be 

3=7 mi m' 

^± = 96m/ dh ± dh/dOm , (67) 
TT^/iOm - l/eOm) = <t/(^- + V^^rr.) 

= e/Oh T Oh . (68) 
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With these new fields and magnetizations, the free energy per site reduces to 

fh = -I {ehOn. + i/ehem) {mh + m)h^) 

+ \ {OhOm - 1/OhOm) {mh) + rh^h). (69) 

As we will show, the choice of the coefficients 6, 9m and 9h will be dictated by physical 
arguments in order to ensure compact and familiar expressions for the critical behavior. 

IV. BINARY LATTICE GASES WITH SHORT-RANGE ATTRACTIVE INTER- 
ACTIONS 

To obtain explicit results for critical behavior we focus now on binary systems with short- 
range interactions (in addition to the hard cores already accounted for). Accordingly, we 
suppose that the small-k expansions of the interactions in k space are 

Zrik) = ZriO) [1 - k^Rl + Oik^)] , (a,r = 1,2), (70) 

with fixed range parameters -R^r- Moreover, to ensure simple criticality, we suppose that 
the interactions are "overall attractive", which we take to mean that jo, AJ12 and A J are 
real and satisfy 

Jo = lJi2(0)>0 and A|Ji2(k)| > 0, Aj(k) > Vk 7^ 0. (71) 

These conditions are easily fulfilled, as, for example, when ^12(1") = -^12 (—r) while Jii(r) and 
J22(r) are positive for all r. 

A. Critical loci 

To identify the singularities of the binary systems, we recall that they are signaled by 
the vanishing of one of the eigenvalues of A, which occurs first when A_(k; A) vanishes. As 
shown in Appendix [Al these singularities arise only when (i) k = 0, and, thence, (ii) when 
the spherical fields A satisfy 

«(Ai,A2) = 0, (72) 

provided the asymmetries of the interactions are not too extreme in the sense that, as we 
suppose henceforth, AJ'''(k) [defined in (I35b[) ] satisfies the conditions (jA2[) and (1A7I1 . 
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Now, any state of the system is specified ab initio by the three thermodynamic fields 
(T, yUi, /i2) which, via (l2Tl) . give A(T, h), and then the densities m. However, for the location 
of critical points, it is more convenient to utilize the set of variables (T, m,m^) introduced 
in ( l63ll and then to solve for (A, A^) as defined in (l59ll . 

Next, let us choose 9m>0 in (i63l ) so that rhc vanishes at criticality or, in other words, 
take 9rn = {—'mi,c/fn2,cY^'^ ■ This condition will be analyzed below in seeking a critical point, 
at a given value of rfiK Likewise, we choose 6* > in (l59l) so that Aj = 0. This condition 
then enforces the link between 9 and m'l", since it implies 9"^ = \i,c/^2,c- As established in 
Appendix [Al the singularities are characterized by u = 0, which via (l62l ) means Ac = jo and 
thence, Ai^c = 9 jo and A2,c = jo/9. 

At this point, one must pay attention to the behavior of integral expressions (l55l l and ( l56l l 
when A_(0) approaches zero. If we accept ([70l) we find that A_ varies as k"^ when A_(0)=0 
and then Q(\) and CaW remain finite at the singularity provided d>d^ = 2 (in this case) 



as seen m 



22|. Accordingly, from here on we suppose the dimensionality exceeds d = 2 and 



may then write 

g{K) = go{9)/jl C,{K) = go{9)LA0)/jo, {(r = 1,2). (73) 

The residual ^-dependence arises from A J and A,P , see f l60l ). (IGTI) , ( l55i l and f l56ll . 

Putting these considerations together we find that the critical locus, Tc(?ia^) with mc = 0, 
may be defined parametrically via 

^"^^^"^"^o(^)[l/^ + /2,o(^)]' ^^^^ 

- ^ (g-l/g)+/i,o(g)-/2,o(g) 

«-l/^^^)+^l,0(^^)^^-/2,0(^)/^^^' ^ ^ 

In fact, the latter relation must be seen as an implicit equation giving ^ as a function of 
while, as shown below, 9m will also be related to 9. Hence, if one realizes that tji^ (i.e., some 
combination of the densities other than the total density) characterizes the composition of 
the system, the function Tc(m^) describes naturally the composition dependence of criticality 
in the binary fiuids. 
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B. Critical neighborhood 



We seek an expression for the physical properties of the binary system in terms of 
(T, m, m^) near the critical locus (Tc, 0, fh^). For this purpose, we first solve for in terms of 
T, rh and fh^: this can be done implicitly in the general case by invoking f l64l) . and explicitly 
in the vicinity of a critical point by implementing a perturbation scheme at fixed rh'' . To 
this end, consider the critical point at Tc(m^) and mc = and its vicinity defined by the two 
small parameters, t a (T — T^), as introduced in (l57ll and m. By construction, u and A^^ 
are small parameters near criticality, so that the integral involved in Q (A) can be computed 



as usual in spherical models, see \4, lO, l22||. However, a significant new feature is that the 
integral is now a function of two vanishing parameters, u and A^. The appropriate extension 
of the standard critical expansion {2^ yields 

GiX) = goJo' [1 - p{l - p^X^y/^ + qou + giX^/jo + 92X^' /fo 

+ 0(^M^/^A^',A^') +o(u)] , (76) 

with coefficients p, p\ Qq, gi, and g2 which in general still depend on 9, and with the critical 
exponent 

7 = max[2/(d-2);l] . (77) 
The integrals Ca{X) are less singular and one finds 

C^iX) = goJo' [L,o + L,iX^/jo + L/X^'/Jo + 0{u, A^')] . (78) 

We note that gi and /^^i vanish in the symmetric case when AJ^ = 0, while gi, g2, p\ la,o, 
la,i and /(j,2 are all of order AJ^^/Jq. 

Using these expansions one can explicitly expand the terms in the spherical constraints 
(l64ll about their values at criticality, which then provides the required relation between u^^'^, 
X^ and t and m. To proceed further, we aim to choose the mixing parameters and, explicitly, 
9m J to ensure that the resulting expansion for u^^'^ begins at orders t and m^, as in (1581 ). 
rather than with m as ( l64ll naively implies. This can be done by imposing the condition 

M ^ l-g,{9)-9h49) 
9[9 + 9g,i9) + h,,{''''' 



To study this, consider first the symmetric case when gi = 1^^1 = 0; the condition then 
reduces simply to 6'^ = 1/6', which, in combination with (l75ll . leads to the equation 

(l + mt'/,,o)(l-^')=0. (80) 
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The conditions fl7T|l enforce AJ^-cr > in symmetric systems which leads to /i,o = ^2,0 > 0. 
Consequently, the only positive solution of this equation is = 1, independently of Thus 
we obtain 6 = 6^ = 1- Finally, we discover, as naturally expected, that symmetric criticality 
is confined to the manifold r/ic = or (pi + p2)c = 1^(7^? while the critical locus is given 
explicitly by the simple parabolic form 



— = l-m' =1-Vq{pi-P2) 

'- c.max 



rpsym ^ ^ "UVKi Fzy 5 (81) 

r,.ma> 



where A;,T|SiJjo = lM(l)[l + /.,o(l)]. 

In the general, nonsymmetric case the condition ([79l) is less tractable and might even 
lead, one could suspect, to multiple solutions. To keep the analysis at the simplest level, we 
note that the coefficients gi{6) and la,i{6) are actually of order Jii(0)/jo and J22(0)/jo- For 
the present work, we will, thus, restrict attention to systems in which the (1,1) and (2,2) 
interactions are sufficiently small relative to the (1,2) attractions (which, then, predomi- 
nantly drive phase separation and yield criticality). In these circumstances the right hand 
side of (TfOll remains positive, ensuring solutions for real 9m and 9: we then select a positive 
root for 9m- 

Supposing then, that the condition f l79l ) is satisfied, the spherical constraint in the general 
case has the expansion 

pu^/^ [l + O (m^"^/^, t, m)] = ctt + Cmm^ + C(tm, m^) , (82) 

where, with the coefficients p, /i,o, ■ ■ ■ defined via the expansions (1761 ) and ( TTSl ), we find 

^* = '+ 99^^ + 1/9 

and, with 

co=jo/9oknTc and w{9) = 1 - gii9) - 91^^. (84) 



2ci 







001 + ^/ooi 

29^9tcl , .2 
w^[9) 



1 + 9« ^^^^ ~ ^ + 9n + o Ol,2 + ^2,2 



(85) 



These expressions are derived only for 7> 1, but the general expansion ( 1821) remains valid 
when 7=1 with, however, different coefficients. For AJ^cr small enough, q and Cm are 
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positive (which we suppose from here on). Thus the structure of (l82ll leads to the usual 
form of the critical singularity in the spherical model. 

The second spherical field A^, which, recalling f l38ll enters into u, is given by 

AVjo = 2eel,fh7h^Co/w{e) + c'^t + d^m^ + mt, m^) , (86) 

where, in similar fashion, the coefficients are found to be 

4 =^2,0 - h,o/9]/wi9)il + l/9^9i) , (87a) 



c. 



2 

/ ^^m 

m 



w{e){i + i/e^e^„ 

1 - l/9'9t + mlco{9h,2 - h,2/e - 2g,)m^' /w^e)] . (87b) 

It should be noted that in the symmetric case, where 6 = 6m = ^, both these coefficients, 
4 and c^, vanish. At this stage, having obtained expansions for u and — and thus for 
Ai and A2 — as functions of T, m, and m^, we are in a position to derive all the physical 
properties of the system in terms of the ffuid variables T, ni, H2, and pi and p2- 



C. Equation of state 

To calculate the equation of state, we need to rewrite the relation (l66l) for the field h using 
the expansions of u and A"^ in terms of t and m, at fixed m\ together with the expansion 

A - Jo = (n + At')/2jo + OiX^\uX^\ u^), (88) 

which follows from (l39l) . At this point, we resolve the freedom to choose 9h in (l65l) by 
requiring that the resulting expression for h is minimally singular. We achieve this by 
canceling the 0{u) term introduced by the factor m^(A — jo) on the right hand side of (1661 ). 
by choosing 

0h = eOm , (89) 
so that u;_ in (l67l) vanishes identically. With this choice the equation of state can be written 

h- \hc + joCht + n7h + 0{mt, f)\ = 

jQ^p~^m{ctt + Cmrh^ + 0{m^,tm)y [l + C>(t, m, m^"^/^)] , (90) 
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provided the expression in braces, which derives from m, remains nonnegative; otherwise, 
this expression must be replaced by zero. Recall indeed, that u must be nonnegative for the 
free energy to be well defined. 

On the left hand of fl90ll the coefficient 



2 

m ' 



(91) 



serves to specify the critical fields, hi^c and /i2,c [via ( l65i l and ( l88l l] and thence the critical 
chemical potentials fii^c and fi2,c- Note that he vanishes with so that in a symmetric 
system, where 6 = 6m = 6h = ^, criticality occurs, as natural, when hc = 0. 
The linear term in t, with mixing coefficient 



Ch{m^) = 2m"'" 4, 



(92) 



similarly determines the near-critical T-dependence of the chemical potentials, /is,i(7';m^) 
and /is, 2(^5 m^), on the phase boundary, a feature to be anticipated in binary fiuid mixtures. 

Finally, note the linear term in m on the left hand side of the equation of state (iQOll : this 
is quite unanticipated from the perspective of previously studied spherical models, at least 
to the authors' knowledge. The corresponding coefficient, which for reasons to be explained 
below, we call the "demagnetization factor," is given by 



(93) 



in which further powers of fh should be noticed. In a symmetric model, this factor simplifies 
to the fairly explicit expression 

+2 



u(m,m1") = 4jo(l + la^ 



m 



1 — m"!" 



2m 



2\2 



(1 - mr) 



+ 0(m^ 



(94) 



Finally, as regards the second external field h\ or chemical potential, near criticality we 
have 



- [hl + n'rh + 0{t^,mt)] 



Jo V '^m^ 



[I + 0{t,m,u 



(95) 
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FIG. 2: (Color online) Phase boundary in terms of the magnetization fh = |(mi + 7712) and m) = 
i(mi — 777,2) in the symmetric case. Criticality occurs on the bold line. The surface represents the 
limit of the single-phase region. Below the surface the parameter u sticks at zero. 

where the right hand side has the same form and is subject to the same conditions as in 
( l90l) while the critical point value is determined by 

hl{rh^)=jom\i + eelYieel, (96) 

and the modified demagnetization factor is 

q2qA _ 

n'(m,mt) =jo^^^ 

+ 4joCo— ^mmt 1 + cq— + O(m^). (97) 
For symmetric models, these three relations reduce simply to 

= 4:]ow) + 0{m,t^). (98) 



D. Phase diagram and critical behavior 

Our result ( l90ll and subsequent relations (l9Tl)-(l93ll. describe the equation of state, i.e. the 
relations between the densities, mi, m2, or pi and p2, and the fields, /ii, and /i2, or chemical 
potentials, /ii and /i2, at temperature T close to criticality. To reveal specific, characteristic 
features we consider, first, the phase boundary in terms of the sum and difference densities m 
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and m'''. It follows from (iQOll that the phase boundary below and up to Tc{m'^) is determined 
by the relation u{t,m,m)) = 0. Figure 2 depicts the boundary in the space (T, m,m^) for 
a symmetric system, for which the critical locus was already derived in (ISTl ). Evidently, at 
fixed m) and for T < Tc(m^) there is a composition gap ^p{T) = Pa — Pi3, where a and f3 
label the two phases. This vanishes as T Tc(m^) and from the magnetic perspective is 
most readily expressed in terms of the spontaneous magnetization which is described by 

mo(T) ^ B\tf with /? = |, (99) 

where i? = (ct/cm)^^^- In fact the critical exponent /9=| represents the standard "universal" 
spherical model result! 

As regards the other critical exponents of the general binary fiuid model, our choice of 
the mixing parameters 9, 9^ and 9^ at fixed ensures that they basically match those of 
the corresponding single-component spherical models. This, of course, is in agreement with 
general considerations of the thermodynamics of multi-component fiuids 3|- Thus regarding 
the density correlation functions, the decomposition (1491 ) shows that the dominant behavior 
of the density structure function near criticality is given by 

^^^(k; T, p) ^ l/{u + k^Rl + ...), (100) 

where Rc is a nonzero range parameter, while on the critical isochore, m = Thc = 0, we have 
u ^ V as follows from (l82l) . Hence, we find S'jvjv(k = 0; T, p^) which is consistent with 

the definition of the critical exponent 7 via, say, light scattering experiments. At criticality, 
this result also implies S'ivjv(k; T^, pc) ~ 1/^^, which confirms the classical value ri = of the 
critical point decay exponent. 

Moreover, in leading order close to but above criticality one can establish the scaling form 

S^^{k;T,p,) ^ t-^X^4k^^{T)], (101) 

where the density correlation length is C,N(T)^^Q/t'^ with critical exponent ^^ = ^7 in accord 
with the general scaling relation 7= (2 — r/)//. The scaling function X^Nix) has the standard 
Ornstein-Zernike form a 1/(1 + x^). 

However, the equation of state in terms of T, m, and the field h requires further exami- 
nation. Thus, while the standard spherical model singularities embodied in (l58ll are evident, 
new features arise from the "demagnetization" term nm on the left hand side of (j90|) . To 
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FIG. 3: Schematic depiction of the equation of state at fixed m' (7^ 0) in a symmetric system for 
temperatures above, at, and below criticality. In the latter case, the equation of state reduces to 
the linear demagnetization form h = nfh (dashed line) for —tjiq < fh < rfiQ (where mo ~ tl^) and h 
in the interval (— /iQ) ^o)- 

understand their significance, consider symmetric systems (i.e., with ^11 = ^22), where (l90l) 
reduces to 

hK.nfh + joP~^{ctt + Cmfh^V, (102) 

while the demagnetization factor n(m,m^) is then given by (i94ll and vanishes only if m^ = 0. 
More generally, however, we note that n{fri = Q,7n^) vanishes in asymmetric systems, only 
on the special two loci fh) = ±mj where one finds 

rhf = (1 - 00l)MO)/^O'Otco. (103) 

It now follows from (|102p that the inverse thermodynamic susceptibility or partial com- 
pressibility, 1/x = {dh/drh)t, does not vanish at the critical point t = 0, Thc = 0: see Figure 
[3l Consequently, the susceptibility does not diverge near criticality in the general case! The 
culprit is clearly the demagnetization effect, i.e., the term nrh which arises both from the 
compositional asymmetry (when fh^ 7^ 0) reflecting the multicomponent nature of the bi- 
nary fluid, and, from the underlying vectorial character of the order parameter: recall that 
spherical models correspond to the n ^ 00 limit of systems of vector-valued spins. Indeed, 
as we demonstrate in the next section, the origin of the demagnetization effect in spherical 
models can be understood directly in terms of vector spin models. 

The nondivergence of the susceptibility/compressibility x means that standard isother- 
mal plots of the chemical potential (or, similarly, the pressure) vs. density or of magnetic 
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field vs. magnetization near criticality take the form illustrated in Fig. [3] with, in general, a 
nonzero slope at T = Tc fixed by the value, ric, of the demagnetization coefficient. Note, in 
particular, that below Tc what would in a standard fiuid system be a constant isotherm of 
zero slope through the two-phase region — i.e., the interval set by the spontaneous magne- 
tization, mo(T) — is now replaced by a straight line with the same fixed slope Uc (at least 
close to Tc). 

It is this fact that leads us to call this anomalous behavior, certainly unphysical in a 
fluid model, a "demagnetization effect". To be more specific, in a real magnetic system with 
long-range dipole-dipole interactions, one must distinguish between the externally applied 
magnetic field -ffext, analogous here to h (or h), and the internal field, ifint, which is what is 
"seen" by individual atomic and molecular spins. The relation between these may be written 

H-.ni = i^ext - iVM, (104) 

where M, analogous here to m (or tti) is the magnetization while is the demagnetization 
2^- More generally the fields Hint and Hext are real-space vectors, as is 



coefficient 
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M, and the relation (I104p can be used only when the system is in the form of an ellipsoid 
and directions parallel to the major axes are considered. (In the case of a sphere one has 



N = ATr/3 [24|,l25|,l26j.) 

One might, in light of these considerations, ask if one should not, similarly, be able to 
introduce an "effective internal field", hi^t ~ /^ext — nm, that would play a natural thermody- 
namic role. However, on the one hand, given the implicit variation of n(?7i, m^^), this seems 
unlikely to be related to the basic thermodynamic parameters, T, /ig-, and per, sufficiently 
directly to be of real value, and, on the other hand, the higher order terms in ( l90l ) indi- 
cate that the linear slope shown in Fig. [3] for the "two-phase" region, will become nonlinear 
outside the critical region. 

Finally, as a further caution, another unphysical feature of the present multicomponent 
luid models must be noted. Indeed, it enters even in single-component spherical models 
! Specifically, whenever 7 > 1 and tie = so that the susceptibility x diverges to 
00 on approach to Tc along the critical isochore, it also diverges when the phase boundary 
is approached below Tc. Even for Uc 7^ 0, a corresponding anomalous feature arises and is 
embodied in Fig. [3] where the slope of the isotherm below Tc remains continuous through 
the phase boundaries (marked by open circles); but in realistic fiuid models there would be 
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FIG. 4: Vectorial representation of the two-species spherical model within a mean-field picture in 
which mi and 012 are the two sublattice magnetization vectors that can rotate with respect to the 
direction, the y axis, set by the parallel external fields, hi and h2. 

breaks in the slope! 



V. VECTOR SPIN MODEL ANALYSIS 



precisely to the n^oo limit of appropriate systems of n-component spins 
character of the order parameter is thus a trademark of the model {4 



To understand the origin of the demagnetization-like effects that enter the present multi- 
component spherical models, it is helpful to recognize that spherical models correspond 

]. The vectorial 
], coupled here 

to asymmetry and multicomponent features. The 'secret' of the demagnetization effect 
appearing first in (l90ll can then be understood by regarding our binary-fluid spherical models 
as magnetic models vv^ith two classes of spins on separate sublattices, just as in ffTll- (fT4ll . but 
nov^ as flxed length vector spins, So^(R°^) (cr = l, 2), rather than scalar Ising spins as originally 
contemplated. 

To obtain insight into the behavior of the model below T^, we may use a simple mean- 
field approach by representing the overall sublattice magnetizations by two mean values, mi 
and m2. The lengths of these magnetization vectors should ideally be taken as m°(T), the 
spontaneous magnetizations (at fixed T<Tc); but it suffices here to consider the symmetric 
situation and so accept equal fixed lengths |mi| = |m2| = 1. 

Then, in contrast to most realistic magnetic systems, it is imperative to allow for the 
imposition of two distinct external magnetic fields hi and h2, corresponding, as fundamental 
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for fluid models, to two distinct chemical potentials /ii and /i2. However, for the "chemical 
interpretation", we must take hi parallel to h2 and may identify the preferred direction as 
the y axis: see Fig. [H The components (mi)^ and (m2)y then correspond to the densities 
mi and m2 in our previous analysis, while 

h = \{hi + h2) with K = {K)y (105) 

describes the external field/chemical potential. On the other hand, h'^ = |(/ii — /i2) char- 
acterizes the compositional asymmetry of the system (even when the two species, here the 
magnetic sublattices, are symmetrically related): that asymmetry is at the heart of the 
matter. As regards the vector-spin dimensionality, however, it suffices to allow for only one 
more dimension and so, regard the as XY or 0(2) order parameters. 

Finally, beyond the symmetric mtrasublattice ferromagnetic couplings (that lead to the 
spontaneous magnetizations), we allow for the intersublattice interactions by a coupling 
constant j > (analogous to jo above). Thus we take the essential part of the mean-field 
free energy to be 

J-'{h; /i^; mi, 1112) = — hi-nii — h2-m2 — jmi-m2. (106) 

Here the external field h is the control variable while is fixed and, as usual, is to be 
minimized with respect to mi and m2. 

Let us, as in Fig. [H introduce the mean tilt angle between the y-axis and m=|(mi + 
m2) and the splitting or separation angle 5, between the m^ and the m. For simplicity we 
suppose that 6 is small (which requires \h^ ^ j); then minimization on 6 yields 

2/7t^siTi^0 9 

^min(e) = -J - 2hcosQ- -—r ^ + 0{h^ If). (107) 

2j + h cos B 

Consider, first, the fully symmetric case in which /i^ = (and 6^i^ = 0). Minimizing 
this expression on 9 then gives Omin = for /i = but Omin = for h <0. This evidently 
corresponds to the usual ferromagnetic situation in which {neglecting dipolar interactions 
and demagnetization effects) the magnetization m switches abruptly from m = {m^)y = —1 
to +1 as the field h passes through zero. 

On the other hand when /i^ 7^ the minimizing value of assumes a nontrivial value for 
h between the limits ±ho given by 



jlJl + h^''/f-i]^lh^'/j^ (108) 
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up to corrections of relative order {hy jY . As a consequence, the magnetization fh(h) no 
longer jumps discontinuously at /;. = from m = — ltom = +l but rather varies continuously 
and almost linearly over the interval —Jiq <h<hQ according to 

m{h) ^ ^ [1 - Wljf] I [I - lih^jfih/hof] . (109) 

More explicitly to leading order in {h^jY one finds for \h\ kJiq, 

2j 



m = cos Gmm , 

n 



(110) 



while for \h\>ho one has cos 9 = sgn{/i}. 

In words, for an external field h not too large compared to the square of the asymmetric 
field , the spins cant themselves in a direction 0min 7^ 0, with, indeed, 9min = vr/2 when 
h = 0\ This minimizing behavior of vector spins is clearly the origin of the seeming demag- 
netization effect and explains our result for the spherical model. Indeed, near the origin, 

2 

for h/ho 0, we find cosGmin ~ 2hj/h'^ which leads to a nondivergent susceptibility of 
magnitude 

X={dm/dh)^-^^^ = 2j/h^\ (111) 

Thus, as the spherical model itself, the asymmetry of the spin model, coupled to the vec- 
torial character of the order parameter, leads to a non- vanishing inverse susceptibility near 
criticality in the general case. Note also that, as in the spherical model, the divergence of 
X re-emerges in the symmetric case when h'^ = 0. Finally, the nonlinear terms in h implied 
1091) s how that one cannot hope to find a simple demagnetization description as in 
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VI. CONCLUSIONS 



We have introduced multicomponent generalizations of the standard spherical model 
that embody lattice-gas hard cores for many-species fiuids by using interlaced sublattices. 
Taking into account a spherical constraint for each distinct species, we have obtained exact 
expressions for the free energy and pair correlation functions in the general case, in terms of 
the basic Fourier space interaction matrix. We have then focused on binary fiuids where the 
diagonalization of 2 x 2 matrices leads to relatively simple results for the physical properties 
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of the system. We find that density and (for ionic fiuids) charge correlation functions can 
be decomposed naturally in terms of two eigenmodes. This formulation, which could well 
lave broader validity, has dramatic consequence for charged fluids as we expound elsewhere 



fly, 
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The present article considers fluids where, in addition to hard cores, only short-range 
overall attractive interactions are present. We show that with an appropriate choice of 
variables (in the form of linear combinations of the mean magnetizations/densities or ex- 
ternal fields/chemical potentials for the two species), the usual critical properties of single- 
component spherical models can be uncovered in accord with general thermodynamic argu- 
ments for multicomponent solutions. Specifically, as the relative composition of the system 
varies, criticality is realized on a well defined locus in the full phase diagram. 

However, an unexpected and intrinsically unphysical "pseudodemagnetization phe- 
nomenon" arises that, except on certain submanifolds, prevents the usual divergence of 
the thermodynamic susceptibilities/compressibilities at criticality. This feature, undesirable 
for model fiuids, is found to be a consequence of an interplay between species and com- 
positional asymmetry and the multidimensional characteristics represented by the hidden 
vectorial character of the order parameter in spherical models. The behavior can, indeed, be 
understood via a simple mean-field description of a corresponding XY spin model. Despite 
this artificial aspect, requiring caution in interpreting results, further aspects of the multi- 
component spherical models seem worth pursuing: on the one hand, ternary fiuid models 
with, say, positive and negative ions in solutions of neutral molecules, could be instructive; 
on the other hand, magnetic systems with different types of ions could reveal interesting 
behavior. 



Acknowledgments 

In the period 2002-2004, this work was supported through the National Science Fondation 
under Grants No 99-81772 and 03-01101. J.-N. A. appreciates both support and hospitality 
from the Institute for Physical Sciences and Technology at the University of Maryland. 



29 



APPENDIX A: LOCATION OF SINGULARITIES 



In this appendix, we locate the singularities of the binary spherical model for the generic 
case of short-range interactions. The singularities derive from the vanishing of A_(k; A) and 
we obtain sufficient conditions to ensure that they arise only when k = 0. 



1. Vicinity of the origin 

We first consider the small-k behavior of A„(k; A). Owing to the second condition (FTTl ). 
the range Ru defined in (jTOl) satisfies Rl2>0. Note also the relations AJ(k) ^ jo-R^^^ and 
AJ^{k)^joR^ k\ where R and R^ are characteristic ranges. At leading order in k we then 
have 

A_(k;A)-A_(0;A) , joR^-X^R^' 



+ + (Al) 

Jo 



which will be valid in a domain |k| < km>0. From the third condition (ITTI ) we find R^>0 
and if, recalling the definition ( 135bll of AJ^, we accept the further condition 

\R^^\ = lim Ajt(k)/jofc' < R^ (A2) 
we are assured that k = is indeed the minimum of A-fk; A) when k < km- 



2. Remainder of the Brillouin zone 

Consider now the subdomain B' of the Brillouin zone B, consisting of all vectors k outside 
the origin domain |k| < km- In the symmetric case when AJ^ = (marked by superscripts 
sym), the second condition fl7T|l leads to |Ji2(k)|<jo so that 

AA!_y™ = A!f"(k; A) - A!_y™(0; A) > A J(k) . (A3) 

But, according to the last member of (j7T]l . there exists a (5A!^™ such that 

AA!_y'"(k; A) > SA'I"^ > , for k G i5' . (A4) 

In the general case where AJ^^ is arbitrary, we may write 

A_(k; A) - A_(0; A) = AA'I'^ik; A) + 5A_[Ajt(k)] , (A5) 
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where from (l36ll we define 

5A_[Ajt(k)] = ^At^ + \JU\.) - ^[At + AJt(k)]2 + 1 J2^(k) . (A6) 

Then, noting that |9((5A„)/9AJt(k)| < 1 and AJt(O) = 0, we see that |5A_[AJt(k)]| < 
|AJ^(k)|. Hence, accepting the further condition 

max[|Ajt(k)|] <5A1^-, (A7) 

which means that the asymmetry is not too strong, one concludes that A_(k; A) > A_(0; A) 
for all k in 
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